Condensate wave function of neutral Bose atoms in power-law traps 
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The macroscopic quantum states of a condensed neutral Bose gas in one-dimensional power-law 
traps are obtained by solving the Gross-Pitaevskii equation numerically. A suitable candidate for 
a trial wave function for the variational calculation of ground state energy is suggested, and shown 
that the suggested function produces a lower ground state energy than the conventional Gaussian 
trial wave function. 
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Bosc-Einstein condensation (BEC) is a first order 
phase transition of bosonic particles when the thermal 
de Broglie wavelength exceeds the mean spacing between 
particles. The first BEC was observed in a series of the 
landmark experiments in 1995 on alkali vapors in which 
the atoms were confined in magnetic traps and cooled 
down to sub /i-Kelvin temperatures. Some of the rele- 
vant features of these trapped Bose gases are that they 
are inhomogeneous, and a harmonic external potential is 
typically used to describe such a trapped system. How- 
ever, due to the development of the experimental tech- 
niques of magnetic traps, the shape of the external trap 
can be managed intentionally, and power-law traps in the 
form of Vtrap ~ r k (k > 0) are accessible now. 

Bosonic atoms in optical lattices and periodic poten- 
tials often require non-harmonic traps such as elliptic or 
a trough type traps 0,0- Therefore, it is necessary to 
predict the condensate wave functions in arbitray power 
traps since it has related with theoretical study of BEC 
such as the variational calculation of the ground state en- 
ergy of the condensed atoms. In harmonic traps a Gaus- 
sian trial function of exp(—r 2 ) is generally applied to the 
variational calculation, but the form is not guaranteed 
the adequate low energy in the power-law traps. Surpris- 
ingly, the shape of the condensate function at power-law 
traps has not been investigated in spite of many stud- 
ies of power-law traps |3J, ky. Most of the previous works 
reported mainly transition temperatures and condensate 
fractions. 

At zero temperature, the condensate wave function 
of the typical model can be described by a non-linear 
Schrodinger equation, named as the Gross-Pitaevskii 
equation (GPE). It is a mean-field approximation for the 
macroscopic wave function of weakly interacting bosons. 



In this short note we'll solve the one-dimensional GPE 
at a few power-law traps numerically, and suggest the 
shapes of the new condensate wave functions that lowers 
the ground state energy than the Gaussian form. Then, 
we will follow the three steps to compare the three en- 
ergies: 1) Apply the above numerical wave functions to 
the Hamiltonian directly and obtain the ground state en- 
ergy numerically. 2) Apply the new trial function to the 
Hamiltonian, and obtain the lowest state energy analyt- 
ically. 3) Apply the conventional Gaussian trial function 
to the hamiltonian and obtain the lowest state energy 
analytically. 

The GPE for N interacting bosons confined by an ex- 
ternal trap potential Vtrap an d under delta-type interac- 
tion is written as 



ih^(r, t) = |-|^V 2 + V trap (r) + g\*(r, i)| 2 J *(r, 

where g is the interaction strength given by g = 47rfi. 2 a/m 
|(| . The a is the s-wave scattering length and m is the 
mass. The *fr(r,i) is real and normalized to the total 
number of particles, N. We choose the power-law po- 
tential for the external trap of the form Vt ra p( r ) = crk 
where c is a positive constant that satisfies the dimen- 
sionality condition, and k = 1,2,3, .... Note that k = 1 
corresponds to a linear trap, and k = 4 to a quadratic 

trap, and so on. 

A scaling of variables such as = y/ h/moj for the 
length, and Hco for the energy reduces Eq. Q to a dimen- 
sionless form 0-0 It corresponds to h — m = lu = 1 in 
Eq. 0]). We consider the one dimensional case because 
it is good enough for the comparison. 
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For convenience we chose c = 1/2 with an adjusted di- 
mension of energy. The P is the dimensionless num- 
ber parameter defined by P = AttN |a|/a/, , which is 
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FIG. 1: Density profile of condensed bosons at one- 
dimensional power-law traps, (a) for linear trap, k = 1, (b) 
for harmonic trap, k = 2, and (c) for quadratic trap, k — 4. 
Solid line for P — 0.1, long dashed lines for P = 1, short 
dashed line for p — 10, and dotted lines for P = 100. 



the only variable that differentiate each atoms. Note 
that \a\/a,ho is order of 1CP 3 for alkali atoms, and then, 
P=l corresponds approximately to N — 10 2 . ip(x,t) — 
ip(x)exp(—ifj,t/H) is normalized to unity, and /i is the 
chemical potential. 

We have carried out a numerical procedure for Eq. J5J), 
and have obtained the static solution if>{x) at various 
powers. Fig. 1(a) is for the linear trap (k — 1), and 
Fig. 1(b) is for the well-known harmonic trap (fc = 2), 
and Fig. 1(c) is for the quadratic trap (k = 4). For small 
number of particles there is little dependence of the traps, 
and the wave function is close to the Gaussian. However, 
for large number of particles the wave function does not 
resemble the Gaussian any more, and it depends seriously 
on the power of the traps. 

In the strongly repulsive limit of P ^ 1, the kinetic 
energy term in Eq. © is neglected and then the solution 
becomes ip 2 (x) = (fi — \x k )/P. However, it is not appli- 
cable for most of large realistic P's. Therefore, we need 
to propose an alternative trial function for a large num- 
ber of particles. We suggest it as ip(x) ~ exp(—x k ) at 




FIG. 2: The lowest energy per particles, (a) k = 1, (b) k — 4. 
Solid line from the new trial function, dashed line from the 
Gaussian trial function, and cross with dotted line from the 
numerical calculation. 



the fc-th power-law traps where k = 1, 2, 3, .... Then, the 
Gaussian trial function is a special case for the harmonic 
traps. 

To check the effectiveness of the new trial function, 
we need to compare the ground state energy densities 
produced by the new ones with two others: One from 
the direct numerical calculation by the Fig. 1, and the 
other from the analytic calculation by the conventional 
Gaussian trial functions. Similar to Eq. @, the ground 
state energy per particle at one-dimension can be written 
in the dimensionless form 0, @ 
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+ x k ip 2 (x)+P^ 4 (x) 



(3) 



One can minimize the energy densities in Eq. J3J using 
the variational ansatz with the new trial function 



ip K (x) = A k e Ki) , 
and with the Gaussian trial function 



(4) 



(5) 



The d is the variational parameter which fixes the width 
of the condensate. is determined from the constraint 
J^ 00 ip 2 (x)dx = 1. Then, A k = ^ k / {4irT (3 / k)d 3 } . 

Substituting ipK,G in Eqs. and JSJ into Eq. iJS), we 
obtain the energies per particle Ek.g/N — /k,g respec- 
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where T is the Gamma function. It is clear that fx 
when k = 2. 

We have determined the minimum of the energy densi- 
ties by solving dfK,G/dd\d c — numerically, and plotted 
fK,G(d c ) as a function of log P. Fig. 2(a) is for the linear 
trap, and Fig. 2(b) is for the quadratic trap. At small 
number of atoms such as logP < 1, the lowest energies 
by the two trial functions show similar values with the 
numerical calculation. On the other hand, for large num- 
ber of atoms the trial function does not always guarantee 



an effective approximation for the ground state. It de- 
pends strongly on the power of the trap. At the quadratic 
trap the two trial functions are still good enough for the 
ground state even for large number of atoms, but at the 
linear trap it is not effective anymore. However, as we 
see, the new trial functions always produce the lower en- 
ergies than the Gaussian. 

We studied the solution of the one-dimensional GPE of 
neutral Bose atoms under power-law traps numerically, 
and suggested a better trial wave function, exp(—r k ), 
beyond Gaussian. We showed that the new trial function 
produce lower energy densities than the Gaussian. It can 
be applicable effectively for condensed wave fucntions at 
any non-harmonic traps. 
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